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^ , Abstract. We consider a model of a quantized fermion field that is based on the 

' Dirac equation in one dimensional space and re-examine how the fermion number 

of the vacuum, or the vacuum charge, varies when an external potential is switched 
on. With this model, fractionization of the vacuum charge has been illustrated in 
the literature by showing that the external potential can change the vacuum charge 
from zero to a fractional number. Charge conservation then appears violated in this 
process. This is because the charge that has been examined in this context is only 
a part of the total charge of the vacuum. The total charge is conserved. It is not 
fractionalized unless the Dirac equation has a zero mode. Two other confusing aspects 
j—i are discussed. One is concerned with the usage of the continuum limit and the other 

with the regularization of the current operator. Implications of these aspects of the 
vacuum problem are explored. 
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^ . 1. introduction 

OO 

The vacuum charge is the fermion number of the vacuum in units of the charge of the 
particle that is associated with the field under consideration. Throughout this paper we 
^ ■ use the words "fermion number" and "charge" interchangeably. Jackiw and Rebbi [1] 
pointed out that, in the presence of a zero mode (a normalizable eigenstate of the Dirac 
equation with energy zero), the vacuum charge of the associated quantum field can 
take half-integral values. Their observation was followed by a large number of papers, 
mostly in the late 1970s and 1980s, in which various such possibilities were explored [2]. 
According to those analyses there are situations such that the vacuum charge can take 
not only half-integral but any fractional, irrational values. The notion of fractionization 
of the vacuum charge has been discussed in many different areas of physics. 

The fractional charges that have been examined so far can be classified into two 
types, A and B, depending on how the charge varies when the external potential (or 
background field) involved is continuously varied. Type A is such that the charge 
remains invariant under the the variation of the external potential. The half-integral 
charge that was found by Jackiw and Rebbi [1] is of this type. There are other examples 
of this type in which the charge can take values ±| and ±| in certain units pj. In 
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type B, the charge (as it appeared in the literature) varies continuously as the external 
potential is varied. In type A it is essential that the external potential has a certain 
topological structure in terms of its large distance behavior. We do not question type A 
which we think is well understood. Rather we focus on type B. In type B it may appear 
as if some topology of the external potential is involved. As we will show in due course, 
however, the topology related to type B is actually a trivial one. 

In this paper we consider a fermion field which is based on the Dirac equation in 
one-dimensional space and re-examine the charge of its vacuum. This model is a typical 
one with which fractionization of the vacuum charge has been illustrated in the literature 
[2]. We assume that there is no zero mode so that the charge fractionization of type A is 
absent. In the fractionization of type B, when an external potential is switched on, the 
vacuum charge changes from zero to a fractional value. It then appears that the charge 
is not conserved in this process. This is because, as we will emphasize, the charge that 
has been examined in this context in the literature is only a part of the total charge. 
Despite the long history of the vacuum problem, its aspect regarding fractionization 
versus conservation of the vacuum charge has not been fully addressed. The purpose of 
this paper is to clarify and explore implications of this and related aspects. 

Throughout this paper we follow Dirac's hole theory. Alternatively one can use 
the second-quantized quantum field theory. We, however, find the hole theory more 
convenient in analysing the structure of the vacuum. The calculation of the vacuum 
charge involves the difference between two divergent quantities. The vacuum charge 
may or may not appear fractional depending on how the difference is calculated. We 
propose a natural way of resolving this ambiguity. This enables us to clearly identify the 
part of the charge that has been left out in earlier analyses. When the external potential 
is switched on the total vacuum charge remains conserved. It is not fractionalized. There 
are two other related aspects that we attempt to clarify. One is concerned with the usage 
of the continuum limit and the other with the regularization of the current operator. 

Before examining the relativistic field with the Dirac equation it would be useful 
to review the time-honored Friedel sum rule for a nonrelativistic Fermi gas with an 
impurity placed in it [I]. This we do in Sec. II with emphasis on charge conservation. 
The impurity produces a potential, V(r), that acts on the particles in the gas. If we 
choose V(r) arbitrarily, it appears as if the charge induced around the impurity can take 
any fractional value. This is related to the usage of the continuum-energy limit. The 
total charge of the gas is conserved. 

In Sec. Ill we examine a model of the Dirac field in one dimension. We see a 
situation very similar to the one found in Sec. II. In Sec. IV we examine the vacuum 
confined in a one-dimensional "bag" of a finite radius. Even such a finite bag is beset 
with an ambiguity that is related to the regularization of the current operator. A 
summary and discussions are given in Sec. V. We suggest an implication regarding the 
charge renormalization in quantum electrodynamics. We use units such that c = 1 and 
h = 1 throughout this paper. 
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2. Friedel sum rule: Nonrelativistic Fermi gas with an impurity 

Consider a nonrelativistic Fermi gas confined in a sphere of a large radius R (in three 
dimensions). We eventually let R — > oo. When an impurity is placed in the gas, it 
induces a change in the charge density distribution. The free gas (without the impurity) 
is described by the Schrodinger equation 

(H o -E)<j>(E,T) = 0, H = p 2 /(2m), (2.1) 

In the presence of the impurity at the origin we have 

(H- E')^(E', r) = 0, H = H + V(r), (2.2) 

where V(r) is the perturbation potential due to the impurity. It is a central, localized 
potential. The (f)(E,r) and ip(E',r) both vanish at the boundary of the sphere, i.e., 
r = R, and are normalized within the sphere. 

The energy levels are discrete. The E and E' are different. If it is attractive 
(repulsive), V(r) shifts the energy levels downward (upward). There is a one-to-one 
correspondence between the perturbed and unperturbed energy levels. Imagine that 
we introduce the perturbation as r)V(r) and let rj increase gradually from to 1. 
Then the n-th E-\eve\ adiabatically goes to the n-th £"-level. We take this one-to- 
one correspondence for granted. 

The density of the gas, defined as the deviation from that of the free gas, is given 

by 

p{r) = £ hK#,r)| 2 - E l^r)| 2 , (2.3) 

E'<E' p E<E F 

where Ep (E' F ) is the free (shifted) Fermi energy. The total charge, which we define as 
the deviation from that of the free gas, is 

Q = 4tt f R r 2 p(r)dr = E " E • ( 2 - 4 ) 

J ° E'<E' p E<E F 

When the perturbation is switched on, the total number of the particles in the gas 
remains the same and hence Q — 0. (See, however, the remarks given towards the end 
of this section.) We divide each of the density and the charge into two parts, 

p{r) = p c {r)+p d {r), Q = Q C + Q d , (2.5) 
Pc(r)= E mE',r)\ 2 ~ E (2-6) 

E'<E F E<E F 

Pd (r) = sign(E' F - E F ) E l^(£>)| 2 , (2.7) 

{E F ,E' F } 



Q c = 4n [ R r 2 p c (r)dr= E " E> ( 2 - 8 ) 

J ° E'<E F E<E F 

Qd = 4tt f R r 2 p d (r)dr = sign(^ - E F ) E ■ ( 2 - 9 ) 

{e f ,e' f } 



Vacuum charge fractionization re-examined 



4 



The radius R is finite, although very large. The energy levels are discrete. The J2ei 
which is the number of energy levels within the specified range, is well-defined and is 
an integer. Hence the Q's defined above are all integers. 

Let the radial part of the normalized partial-wave function in the presence of the 
impurity be ipi{E',r). In the asymptotic region where V(r) is negligible, we obtain 



1 \ 1/2 1 

ibAE', r) = ( - sin 

\2txRJ r 



k'r + m{E') - y 



(2.10) 



where k' 2 /(2m) = E' . In Eq. (I2.6p . the difference between \ip(E' ,r)\ 2 and \(p(E,r)\ 2 
is small at large distances. Hence the distribution p c (r) is mostly concentrated around 
the impurity. By using the Schrodinger equation and Eq. (12.101) and by taking the 
continuum limit, that is, by replacing J2e with 2(2R/ir) + 1) / dk, one can derive 
the well-known formula [I], 

Q c =-J2(21 + 1)Vi(Ef). (2.11) 
n i 

In deriving Eq. (12. lip one does not have to know the wave function in the vicinity of 
the impurity. If we choose V(r) arbitrarily, the above Q c can take any fractional value. 
This contradicts what we stated at the end of the preceding paragraph. We will discuss 
this contradiction shortly. For a real electron gas in a metal, only if V(r) is chosen 
judiciously, Friedel's self-consistency condition or the Friedel sum rule 

Q C = Z (2.12) 

is satisfied. Here Z is the (excess) charge of the impurity and is an integer. 

Next let us turn to Pd{f)- The Fermi momenta k F of the free gas and k' F of the 
perturbed g related by k' F R + i](E' F ) = k F R and 

k F = k' F -k F = -r ]l (E F )/R, (2.13) 

Unlike k F , k' F depends on /. If the impurity is introduced adiabatically, the particle 
population in each of the partial waves remains the same. The difference E F — E' F 
is very small (because R is very large). If we replace J2e °f Eq. (12.71) with 
2(2R/7r)Ak F £j(2Z + 1) / 8(k - k F )dk, we obtain, 

Pd (r) = --J2^ + lME F )^ 2 (E F ,r), (2.14) 
n i 

Q d = --J2( 2l + 1 ME F ). (2.15) 
71 i 

This confirms Q = Q c + Qd = 0. To conserve the total charge of the gas, it is crucial 
to include Pd(r) which stems from the change in the Fermi energy. In the asymptotic 
region where V(r) is negligible, we obtain 

Anr 2 p d (r)= - A £( 2 / + l)^(^) 



x sin 2 



. . I TV 

k F r + r]i{E F ) - — 



(2.16) 
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This distribution is diffuse; it spreads more or less uniformly over the entire asymptotic 
region. As R — ► oo the density pa( r ) vanishes but Qd does not. 

We noted below Eq. (12.91) that Q c and Qd are integers. This is in contradiction 
with Eqs. (12.111) and (12.151) which show that Q c and Qd can be nonintegers and their 
dependence on V(r) is continuous and smooth. The origin of this contradiction can be 
traced as follows. In deriving Eqs. (12.111) and (12.141) . J2e is replaced with (-R/vr) / dk. 
This effectively smears out the discrete energies. The number of such smeared-out 
energies in a specified energy range becomes fuzzy. In Eq. ( \2.9\\ we count the number 
of energy levels in interval AEp = E' F — Ep, which is of the same order of magnitude as 
the spacing of the discrete energy levels. In terms of k, the level spacing is Ak = ir/R 
which is comparable with Akp = —rj/R of Eq. (I2.13I) . In such a case, smoothing the 
energy levels results in a considerable uncertainty. For the total charge Q, we count the 
number of states in interval {0, Ep} which is much wider than the energy level spacing. 
This is why the smoothing does not affect Q (= Q c + Qd) but it obscures each of Q c 
and Qd- The fractional values that appear through Eqs. (I2.11I) and ( 12. 15ft are artificial. 

The model that we have examined above is just a mathematical model. Let us 
examine its relevance to a real electron gas in a metal. For example, the impurity can 
be a Zn ion of charge +2 which replaces a Cu ion of charge +1 in a metal lattice. The 
excess charge is Z — 1 . In a real electron gas we have to remember the electron-electron 
interaction which is repulsive. This pushes the electrons in the diffuse distribution 
towards the surface of the metal. No diffuse distribution remains inside the metal. 
(This is a well-known result of electrostatics.) 

We have assumed that the impurity potential V(r) is switched on adiabatically. 
This, however, is not what actually takes place in the real electron gas. As pointed out 
by Friedel, when one of the ions is replaced by an impurity of excess charge Z, additional 
electron (s) are added to the system at the same time. The charge of the added electrons 
exactly cancels Z and the entire system, including its background ions, remains neutral 
[HE]. The Fermi energy of the gas with the added electron(s) becomes the same as the 
free Fermi energy Ep. If the electron-electron and electron-ion interactions are taken 
into account in a self-consistent manner, e.g., by means of the Hartree-Fock method, 
V(r) will adjust itself such that Q c = Z is satisfied and the impurity charge is screened. 
Since Z is an integer, no fractional charge appears. 

We raised a question regarding the continuum limit of the energy levels. We, 
however, do not propose to dismiss the Q c and Qd obtained in the continuum limit. It 
is conceivable that Eqs. (12. lip and (12.151) are practically acceptable or even preferable. 
The reason is as follows. Imagine that V(r) is varied continuously and one of the E' 
levels crosses Ep. Then Q c of Eq. (12. 6ft suddenly changes by 2(2/ + 1) where I is the 
angular momentum of the E' level. This is probably unrealistic. The energy levels near 
the Fermi energy may be only partially occupied. The fuzziness that is introduced by 
smearing the energies may be appropriate in such a situation. It is well-known that the 
Friedel sum rule ( 12.121) is useful [H [5] . 
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3. Dirac field in one dimension 

We examine the vacuum charge of a Dirac field in one dimension. We start with the 
Dirac equation for a particle in a given external potential in single particle quantum 
mechanics. We compare the free and perturbed systems. The Dirac equation for the 
free system is 

{Hq - E)(j)(E, x) = 0, H = ap + (3m, (3.1) 

where p = —id/dx, a and (3 are the 2x2 Dirac matrices and m is the mass of the 
particle. In explicit calculations we use a = cr 2 and (3 = cr 3 . The Dirac equation for the 
perturbed system is 

(H — E')ifj(E' J x) = 0, 

H = H + U(x), U(x) = j3S{x) + V(x), (3.2) 

where S(x) is a Lorentz scalar and V(x) the zero-th component of a Lorentz vector. 

If V(x) = 0, there is a symmetry between the positive and negative energy spectra: 
i/j(E', x) and ip(—E', x) are related by if)(—E', x) = a/3ip(E', x). Because there cannot be 
energy-degeneracy in one dimension, this symmetry may seem to exclude the possibility 
of a zero mode. A zero mode, however, is allowed provided that S(x) has a topology 
such that m + S(oo) and m + S(—oo) are of opposite signs [H [6]. Such a situation 
arises with a kink soliton background. The zero mode remains no matter how S(x) is 
modified as long as this topology is maintained. If V(x) ^ 0, the symmetry between 
the positive and negative energy spectra is broken. But the zero mode can still exist 
as long as [m + S(oo)][m + S(—oo)] < and |V(:r)| < |S(a;)| when x — ► ±oo. The 
zero mode is stable against the variation of V(r). This is the type of zero mode that 
underlies half-integral charges pQ. 

In this paper we assume that S(x) has no such topology as described above. To be 
more explicit, we assume that [m + 5(oo)][m + S(— oo)] > 0. Under this assumption for 
S(x), it is still possible to have a zero mode. This can be done by carefully adjusting V(x) 
such that one of the eigenvalues becomes zero. This zero mode is obviously unstable 
against the variation of V(x). Its energy eigenvalue can deviate from zero even if the 
variation of V(x) is very slight. We may say that the existence of this zero mode is 
accidental. This is in contrast to the zero mode that we described in the preceding 
paragraph. In any case we assume that our Dirac Hamiltonian H has no zero mode. 

The H of Eq. (13. 2p can be transformed into different forms. An example is, 

, T/ / \ df(x) 

= ap+ V[x) 

dx 

+ P[m + S(x)][coB2f(x) -iasm2f{x)}. (3.3) 

The term with —i(3a = —ia^a2 = —<J\ is a pseudoscalar potential. Potential V(x) can 
be eliminated by choosing f(x) such that df(x)/dx = V(x). If m + S(x) = 0, then the 
transformed Hamiltonian Hf simply becomes Hq = ap, which is of course solvable. We 
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will make use of this transformation in the next section. In this section we assume that 
the mass is nonzero. 

Suppose we eliminate V(x) in Hf by means of df(x)/dx = V(x). It is understood 
that m + S(oo) and m + S(—oo) are of the same sign. Assume that /(oo)/(— oo) < 0. 
Then the peudoscalar part of the potential of Hf, which is proportional to ia/3 sin 2f(x), 
has a topology, that is, its asymptotic values for x — > ±00 are of opposite signs. If we 
go back to H, however, there is no such topology because df(x)/dx has the same sign 
at x — > ±00. Hamiltonians of the form of Hf are extensively used in the literature 
[2] and one might have the impression that the fractional charge associated with such 
Hf is due to its topology. This topology, however, is a trivial one in the sense that it 
can be transformed away. Let us emphasize that the topology carried by the Lorentz 
scalar potential remains the same in transformation ( 13. 31) . The topology of S(x) that 
we mentioned in the paragraph below Eq. ( 13. 2\\ cannot be transformed away. 

We return to Hamiltonian H of Eq. (13.21) . We assume that the system is confined 
within a box of large radius R, i.e., r = \x\ < R. We do this by assuming an infinite 
square-well potential of the Lorentz scalar type, which is 00 for r > R. It is understood 
that this square-well potential is a part of S(x). The wave function has to satisfy a 
boundary condition at r = R; see Eq. (14.11) of the next section. Since R is finite, the 
energy levels are discrete. We let R — > 00 in the end. (We keep R finite in Sec. IV.) 
The wave function vanishes for r > R. Therefore, S(x) and V(x) for r > R do not 
appear in the following calculations. We assume that S(x) and V(x) within the box are 
both localized around the origin. In this way we can see a parallelism between this and 
the preceding sections more easily. In the following calculations, however, the S(x) and 
V(x) within the box can be chosen arbitrarily. 

We can find a one-to-one correspondence between the free energy levels and the 
perturbed ones as we did in Sec. II. For simplicity we assume that positive Es correspond 
to positive E's. In other words, the perturbation does not change the sign of the energy. 
The vacuum is such that the negative energy levels are all occupied. The number of the 
negative energy levels is of course infinite. We, however, first assume that the negative 
energy levels are filled only down to a certain Fermi energy Ep < for the free vacuum 
and similarly to E' F < for the perturbed vacuum. Eventually we let the Fermi energies 
tend to —00. Unlike the Fermi energy of the nonrelativistic Fermi gas, the Fermi energy 
of the vacuum is a mathematical device. We need such a device so that we can clearly 
keep track of the number of particles in the vacuum. 

The density and the charge of the perturbed vacuum are respectively defined as the 
deviations from those of the free vacuum, 



p(x) 



E W£»i 2 



E \<f>(E,x)\ 2 



(3.4) 



0>E'>E' F 



0>E>E F 




0>E'>E' F 



E 



E 



(3.5) 
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The one-to-one correspondence between the free and perturbed energy levels implies 
Q = 0. When the perturbation is switched on, the total charge remains conserved. 

Exactly in the same way as we did in Sec. II we divide each of the density and the 
charge into two parts, 

p(x) = p c (x) + pd(x), Q = Q C + Qd, (3.6) 
Pc (r)= £ mE',x)\ 2 - £ |<K£,*)| 2 , (3.7) 

0>E>>E F 0>E>E F 

Pd (r) = sign(E F - E' F ) E \^E\x)\\ (3.8) 

{E F ,E' F } 

q c = [ R Pc (x)dx = E - E > ( 3 - 9 ) 



" R 0>E>>E F 0>E>E F 

f R 

Qd= Pd(x)dx = sign(£ F - E' F ) E 



(3.10) 



The Q can be rewritten as 



Q = \\ E - E I- (3-ii) 



v 0>E'>E^ 0<E'<^ / 



Imagine the positive-energy level of = \Ep\ of the free system. The E F > of Eq. 
( 13. lip is the energy level (of the perturbed system) that corresponds to \Ep\. Figure 
1 schematically shows how the Fermi energy shifts when the perturbation is attractive. 
In general E' F ^ —E' F . In deriving Eq. ( 13. lip we used 

E + E = 2 E > (3-12) 

0>E'>E' F 0<E><E F 0>E>E F 

which is based on the one-to-one correspondence between Es and E's. Furthermore the 
Q c can be written as 

Qc = l( E - E V (3-i3) 

Z \0>E'>E F 0<E'<\E F \J 

The Q c is a measure of the asymmetry between the positive and negative energy spectra. 



\Ef\ 



E = 



> 



Ep < 7 

1 E'p < 

Fig. 1 Unperturbed and perturbed Fermi energies. 
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When Epi E' F and R are finite, the X)b' s are a U well-defined. The above Q's are 
all integers. When Ep — > — oo, — > — oo and R — > oo, the Q's remain as the same 
integers. In Sec. II and so far in this section we have distinguished the unperturbed 
and perturbed energies, using different notations E and E' . In the following, when the 
distinction between E and E' is unimportant or is clear from the context, we may denote 
both of them as E for notational brevity. 

The p c (x) is the one that has been extensively examined in the literature, for 
example, in Refs. [7J E]. According to those analyses it is given by 

Q c = --v(-oo), (3.14) 

7T 

where we have taken the limit of Ep — > — oo. The phase shift r)(E) pertains to the 
transmission coefficient T(E). In the Born approximation the phase shifts are given by 

r](E) = -\ ( R \EV(x) + mS(x)]dx, (3.15) 
k J~r 

which becomes exact when \E\ — ► oo. Thus we arrive at 

1 r R 

Q c = -- J R V{x)dx. (3.16) 

The S(x) has no effect on Q c . It follows from Eq. (13.151) that 77(00) = — ry(— 00) and Eq. 
( 13. 14)) can be written as Q c = (1/2%) [77(00) — r)(—oo)). This shows that Q c is a measure 
of the asymmetry between the positive and negative energy phase shifts. Let us add 
that, in the nonrelativistic version of the present model, the phase shifts in the Born 
approximation are r](E) = —(m/k) J_ R V(x)dx. Unlike its relativistic counterpart, the 
nonrelativistic i](E) vanishes as E — > 00. 

No attention seems to have been paid to pd(r) in the literature so far. In order 
to see similarity to the calculation of Sec. II clearly, let us assume that U(x) of Eq. 
( 13. 2p is an even function of x so that parity is a good quantum number. There are two 
partial waves, one with even parity and the other with odd parity [9J. The transmission 
coefficient T(E) is related to the partial wave phase shifts f]± by T = [e 2lv+ + e 2ir, -}/2 
and hence 

V (E) = r 1+ (E) + r 1 ^(E). (3.17) 

In the asymptotic region where U (r) is negligible, the perturbed wave function of energy 
\E\ > m and positive parity takes the following form, 



1 (t? \ E + m cos(kr + 7] + ) \ 

where k = ^ E 2 — m 2 and x = x/r. Suffix + of the phase shift rj + (E) refers to positive 
parity. There may be bound states with l^l < m. For negative parity we obtain 



, , \E — m( -E^—xsm(kr + ri^ \ 

ib(E,x) = \ — E ~ m „ v Y J ). 3.19 
r v ' ; V 2RE \ cos(kr + r]_) J v ' 

The 4>±(E,x) are obtained by dropping rj±. 
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The Fermi momenta before and after perturbation U(x) is switched on, kp and k' F , 
are related again by Eq. (I2.13p . By taking the continuum limit, we obtain 

Pd(x) = -J2m(E F M(E F ,x)\ 2 . (3.20) 
n i=± 

By integrating this with respect to x and letting Ep -+ -oowe obtain 

Q d = -v(-oo) = - [ R V(x)dx. (3.21) 

7T 7T J-R 

In the asymptotic region where U(x) is negligible, in the limit of Ep — > — oo, we find 
that \ip + (E,x)\ 2 -> 1/(2/2) and hence p d (x) -> (l/7ri?)r/(-oo). 

The remark given in Sec. II regarding Q c and Qd also applies to the Q c and 
obtained above. As a model, V(x) can be chosen arbitrarily. Then Q c of Eq. (13.161) 
and Qd of Eq. (13.211) can take any fractional values. As we noted below Eq. (13. 13ft . 
however, the Q's of Eqs. (13.91) and (13.101) are integers. The fractional values of Q c of Eq. 
(13.16!) and Q d of Eq. (I3.21j) are mathematical artifacts that stem from replacing discrete 
summation J2e with the integration (R/ir) J dk over the smoothed energy levels. 

At the end of Sec. II we suggested a practical justification for Q c of Eq. (12.111) of 
the nonrelativistic electron gas in a metal. Such justification is untenable for the present 
relativistic model in which the Fermi energy is a purely mathematical device. Towards 
the end of this section we argue that, although they should not be taken literally, these 
artificial results do no harm in real physics. 

Let us briefly discuss an explicit example. Assume that 

S( x ) = 0, V(x) = -2X5(x). (3.22) 

This is a special case (Ai = A2 = A) of the model that was examined in Ref. [10] . This 
model is related to the model examined by MacKenzie and Wilczek [7] by transformation 
(13. 3p . In dealing with the 5-function potential, we start with a square- well potential 
and after solving the Dirac equation we take the narrow-width limit of the square well- 
potential. In this connection see Ref. [TTj . The phase shifts are given by 

tanr] ± (E) = E ^~ m t&n\. (3.23) 
k 

If g < there is a bound state of positive parity. Its energy is 

£ = mcos(2A). (3.24) 

The Dirac matrices used in Ref. [TU] are a = 0\ and (3 = 03, which can be transformed 
to a = 02 and (3 — (73 that we are using by the unitary transformation exp(z7ro"3/4). 
The wave function of the bound state is 



H X)= J«B+Il( 1 W (3.25) 



2m \ -f^x 

■ E+m 



where n = \/m 2 — E 2 . If g > there is a bound state of negative parity but we do not 
delve into such details in this paper. The p c (x) was worked out explicitly in Ref. [10J. 
The summation J2e was done as (R/tt) J dk. The p c (x) is indeed concentrated around 
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the origin. It decays like e~ 2mr . If we let m — > 0, p c (x) becomes uniform. By carrying 
out the integration / p c (x)dx explicitly, one obtains Q c of Eq. (I3.16p . The Pd(x) was 
not considered in Ref. pH] but it is given by Eq. ( 13.201) . It is uniform. Let us add 
that a model in which the (^-function potential is replaced by a square- well potential has 
been also examined [12]. The square well- version is related to the model of Ref. [13J by 
transformation (13. 3p . This ends the discussion of the example. 

As we pointed out already the fractional charges of Eqs. (I3.16P and (13. 2 1 p are 
mathematical artifacts. Still let us speculate on what happens if we accept the fractional 
charges literally. In the model that we have examined, the particles individually interact 
with the external potential but there is no interaction between the particles. We assumed 
that the external potential can be chosen arbitrarily. In real physical systems, however, 
the particles generally interact with each other. When the inter-particle interaction is 
(approximately) eliminated in a self-consistent manner, the potential that emerges is 
not an arbitrary one. If the particle of the field is the electron in a metal, the remarks 
given at the end of Sec. II apply. The repulsion between electrons pushes the electrons 
in the diffuse distribution towards the edge at r = R. Such a charge accumulated at 
the edge can easily move on to another object that comes into contact with the metal. 
For Q c , the Friedel sum rule (I2.12p will be satisfied. Thus the self-consistency prevents 
fractionization of the vacuum charge. 

Our model can be regarded as a one-dimensional simulation of the relativistic 
nuclear shell model. In that case we have a number of nucleons in the positive energy 
levels, but let us focus on the vacuum effect. According to the Dirac phenomenology 
of nuclear physics, the relativistic potential of the nuclear shell model is of the form 
of U(r) = j3S(r) + V(r) [14J. For medium to heavy nuclei, S(r) and V(r) are well 
represented by the Wood-Saxon form. Inside the nucleus S + V ~ —50 MeV and 
— S + V ~ 800 MeV. The potential acting on the lower component of the Dirac wave 
function is — S + V, which is almost as large as the nucleon rest mass of 940 MeV. 
Consequently relativistic effects can be appreciable even at very low energies. If we 
assume V = 400 MeV, we obtain Q c = 400(2R /ir) MeV-fm where #0 ~ (1.2) A 1 / 3 is 
the radius in fm of the nucleus of mass number A. If we assume A = 100, for example, 
we find Ro ~ 5.6 fm and Q c = —Qd ~ —7.2. Approximately 7 nucleons are moved from 
the concentrated region to the diffuse distribution. This number is not negligibly small. 

In our model calculation the diffuse distribution is uniform in the asymptotic region 
(outside the range of the shell model potential). Unlike the electron case, the interaction 
between the nucleons is attractive. When the nucleons are treated in a self-consistent 
manner, like the Hartree-Fock method, the interaction between the nucleons will produce 
an attractive potential and the diffuse distribution will be pulled back to a finite region. 
It may become something like a nuclear halo. The nucleons, including those in the diffuse 
distribution, will form a nucleus of a finite size. The total charge of the nucleus is an 
integer. It is an interesting challenge to do such a self-consistent calculation including 
the vacuum effect. If the nucleons in the diffuse distribution somehow escaped to infinity, 
the nucleus would be left with Q c . If this Q c takes a fractional value, the nuclear charge 
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becomes fractional. This would lead to bizarre consequences. We know that this is not 
the case in reality. 

Finally let us mention super-heavy quasimolecules that could be formed in a 
collision of two very heavy ions. This subject has been extensively discussed by Greiner 
et al. [T5j . When the strength of the Coulomb potential due to the merging nuclei 
exceeds a certain critical value, the lowest positive energy level of the electron dives 
into the negative energy sea. This may lead to the creation of an electron-positron pair. 
Suppose that the vacuum charge is fractionalized in this process. Then can a fraction 
of the pair be created? This does not seem to make sense. 



4. Vacuum confined in a finite bag 

We examine a one-dimensional version of the so-called MIT bag model [16]. We also 
discuss the chiral bag model [TF] towards the end of this section. We keep the radius 
R of the bag finite. In Sec. Ill we divided the density into two parts, concentrated 
and diffuse. Such division is not very meaningful when R is not very large. In Sees. II 
and III, we pointed out that replacing J2e with (R/it) J dk leads to fractional charges. 
When R is finite and not very large, one may think that there is no room for such 
complication and no problem arises regarding charge conservation because the charge 
cannot escape to infinity. Nevertheless, the vacuum charge may be fractionalized in the 
bag, depending on how it is calculated. 

The Hamiltonian of the bag model is the same as that of Sec. III. We assume an 
infinite square well potential for r > R so that the particle is confined to r < R. This 
leads to the boundary condition at r = R, 

(1 - iaf3x)^){x) = (1 + axx)^{x) = 0, (4.1) 

which is equivalent to Eq. (14) of Ref. [18]. This condition makes the scalar density 
ip* 'flip vanish at r = R. 

We begin with the simple case of m = 0. The solutions of the free Dirac equation 
can be classified in terms of parity. For the unperturbed system with H , we obtain 

* +W=iV (-sin fa )< (4 ' 2) 

where suffix ± refers to parity and N = l/^/2R is the normalization factor. The density 
\<f)(E,x)\ 2 = 1/(2R) is a constant. The k and E = ±k are determined by the boundary 
condition (14. IB . The positive energy eigenvalues are given by 
f ,\ (2n — 1W 

n 4R ' n=1 ' 2 ' 3 '---' ( 4 - 4 ) 

where superscript (±) refers to the sign of the energy. All energy levels are equally 
spaced with interval ir/(2R). 
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We now examine the perturbed system with the H = H + V(x). We assume 
S(x) = for simplicity. The H and H can be related by transformation (13.31) . The 
following two wave functions satisfy the perturbed Dirac equation, 



and if)-(x) = —iaip+(x). In order for the above to satisfy the Dirac equation, V(x) does 
not have to be an even function of x. If V(x) is not an even function, parity is not a 
good quantum number. We still retain the suffixes ± because the above solutions are 
respectively related to the free solutions <f)±{x). The above is a generalization of the 
solution of Ref. [19] in which V(x) = \x was assumed. 

The perturbed density is \ip(E',x)\ 2 = 1/{2R). Imposing the condition (14.11) we 
find that the perturbed energies and the energy shift, 



This AE has two remarkable aspects, (i) The AE is of the form of the first-order 
perturbation; yet it is exact. This is so no matter how strong V(x) is. Higher order 
terms are absent. If V(x) is an odd function of x, then AE = 0. (ii) This energy shift 
AE applies to all energy levels, that is, all energy levels are shifted by exactly the same 
amount. 

In the vacuum all negative energy levels are occupied. Start with the free 
vacuum and switch on perturbation V(x). Assume that V(x) is attractive so that 
AE < 0. Recall that all energy levels are equally spaced with separation tc/(2R). If 
\AE\ < tt/ (2R), the lowest positive energy remains positive and hence the number of the 
negative energy levels remains the same. There is no change in the vacuum charge, i.e., 
Q = 0. If 7r/(2R) < \ AE\ < n/R, the lowest positive E becomes negative E'. Therefore 
the number of the negative energy levels increases by one. This new negative energy 
level is empty. If we fill this level by hand, the charge of the new vacuum becomes one 
unit greater than that of the free vacuum. No fractional charge appears. 

Suppose that \AE\ ^> ir/(2R). Then a large number of positive energy levels move 
into negative energy sea. The number of such levels can be estimated as 



We have put the negative sign in the left hand side of Eq. (14. 7p because the number 
should be positive when AE is negative. It is interesting that the above result is 
consistent with the Q c of Eq. (13.161) . However, the number of the levels that goes from 
positive to negative should always be an integer. Equation (J4TTJ) is only a fuzzy estimate 
of the discrete number. 

Equation (13.1 ip for Q also holds in the present case. When E' F — > — oo and 
Ep — ► oo, Q is the difference of two divergent series. The difference depends on how 
the summation is done. We propose to do this in the same way as we did in Sec. III. 
We first keep E' F and E" F finite. Then the one-to-one correspondence between the Es 




(4.5) 




(4.6) 




(4.7) 
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and E's dictates that Q = 0. This remains so in the limit of E' F — > — oo and E' F — > oo. 
This we believe is a most natural prescription for dealing with the series. 

There is another prescription that has been widely used |2Q] . This is to rewrite Q 

as 



where the summation is extended to infinity. (No Fermi energy is assumed.) This 
prescription is based on the so-called split-point regularization of the density operator 
in quantum field theory [21]. Let us apply this method to our bag model. With 
= E n {±) + AE where E n {±) is that of Eq. flUD we obtain 



and similarly for the other summation. By taking the limit of s — > +0 and using Eq. 
(14. 6p we obtain 



We have three remarks to make on the above result. 

(i) Equation (14.101) follows no matter how small R is and hence no matter how large 
the level spacing is. If the level spacing is very large and V(x) is not very strong, all 
positive energy levels remain positive. The number of the negative-energy levels does 
not change. Even then Eq. (14.101) implies that the vacuum charge changes, which we 
find unphysical. The current is zero everywhere in stationary states. The charge cannot 
escape from the bag and hence it must be conserved. 

(ii) Imagine that the energy levels are continuous and replace the sum over n = 1, 2, 3, . . . 
of Eq. (14. 8p with an integral J°° dn. We still obtain Eq. (14.101) . (This result is 
independent of the lower end of the integral.) In this sense the prescription (14.81) is 
equivalent to smearing out the discrete energy levels. It is interesting that this smearing 
leads to the same fractional charge as that of Q c of Eq. (13.161) which is also related to 
the smearing of the energy levels as we discussed in detail in Sec. II. We think this is a 
mathematical artifact. 

(iii) Divide Q into Q c and Qd as we did in Sec. III. (This division is not interesting in 
the sense that, when m = 0, p c and pd are both diffuse.) If we use Eq. ( 14. 8ft . we find 
that Qd disappears. The prescription (14.81) is equivalent to ignoring the change in the 
Fermi energy in our proposed scheme. In the case of Sec. Ill, Q c and Qd may both 
be fractionalized but they cancel each other, Q = Q c + Q d = 0. In the present case, 
fractionalized Q c remains because Qd is missing. The total charge is not conserved. 
This is unphysical. 

So far we have assumed m = 0. We now assume that m > 0. The Dirac equation 
can be solved exactly for Hq but not for H in general. As \E'\ becomes much larger 
than m, however, the effect of the finite mass becomes negligible. For such E's the 




(4.8) 




(4.9) 




(4.10) 
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energy shift is given by Eq. (I4.6p . The one-to-one correspondence between the -E's and 
E's implies Q = 0. Here we are again assuming that positive _E"s correspond to positive 
E's, for simplicity No fractional vacuum charge appears. 

Let us see what happens if we use prescription (14.81) when m > 0. The Q is 
essentially determined by the energy shift for very large E's. This situation is similar 
to the one we saw in Sec. Ill where Q was determined by r](±oo). It is not difficult to 
see by using Eq. (14. 8[) for the cases of m > and m = and remembering that the E' 
becomes independent of m as \E'\ — > oo, that 

Q(m > 0) -Q(m = 0) = 0. (4.11) 

Hence Eq. (14.101) also holds even when m > 0. 

We now turn to the chiral bag model [TTl 1201 [22] ■ Fractionization of the vacuum 
charge has been extensively discussed for the chiral bag model, but never for the MIT 
bag model that we are using. This may give the wrong impression that fractionization is 
a peculiar feature of the chiral bag model. As far as the vacuum charge associated with 
the bag is concerned, however, there is no essential difference between the two types of 
the bag models. 

It is instructive to examine the one- dimensional version of the chiral bag model 
that was proposed by Zahed [22J. The quark is assumed to be massless. The quark 
wave function q(x) for a stationary state obeys the Dirac equation 

(op - E)q(x) = 0, (4.12) 

inside the bag and the boundary condition at x = ±R 

[cos 6{x) — ia sin 6{x) + a^x\q{x) = 0, (4-13) 

where a = a-i and 6{x) is a function that pertains to the soliton with which the quark 
interacts. If we write q(x) as 

q(x) = exp{(i/4)[20(x) - n}a}ip(x), (4.14) 

Eq. (14.121) becomes 



1 dd(x) „ 



ip(x) = 0, (4.15) 



and the boundary condition (14.131) is reduced to Eq. (14. II) . As far as the quarks in the 
bag are concerned, Zahed's model is exactly equivalent to the one-dimensional MIT bag 
model with external potential V(x) = (l/2)d6(x)/dx and the quark mass m — 0. The 
energy spectrum and the energy shift that Zahed obtained are the same as those of of 
the MIT bag model. This situation also holds in three dimensions, that is, the chiral 
bag model can be transformed into the MIT bag model with an external potential. 

From our point of view, there is no anomalous, fractional baryon number associated 
with the bag. We can still have a soliton which interacts with the bag, and the bag plus 
the soliton can be interpreted as a baryon. The baryon number carried by the bag, as 
we see it, is an integer. There is no need to invoke the interpretation that a part of the 
baryon number is carried by the soliton as is done in the chiral bag model [221 123] . 
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We should add that not all the people who worked with the chiral bag model 
assumed that the baryon number in the bag is fractional. In the models for baryons 
developed in Ref. [24J, the baryon number of the chiral bag is unity. Examining the 
vacuum charge of a model similar to that of Zahed, Banerjee argued that "there is 
nothing wrong about departing from the split-point form" in regularizing the current 
operator [25]. We concur with him in this respect. Banerjee then suggested a 
regularization scheme such that, even in the presence of spectrum asymmetry that is 
induced by an external perturbation, one uses e ±sE instead of e ±sE ' for Q of Eq. ( 14. 81) . 
This ad hoc scheme leads to Q = 0. We believe, however, the scheme with the Fermi 
energy for the vacuum that we have introduced in this paper is more natural than what 
Banerjee suggested. 

5. Summary and Discussion 

We have examined the vacuum charge for a fermion field model based on the one- 
dimensional Dirac equation with an external potential. This model is a typical one 
with which fractionization of the vacuum charge has been illustrated in the literature. 
We classified the fractional vacuum charges so far examined in the literature into two 
types, A and B. When the external potential of the model is continuously varied, 
the vacuum charge of type A remains unchanged, whereas that of type B varies 
continuously. We assumed that there is no zero mode for the Dirac equation so that 
the charge fractionization of type A does not appear. We have focussed on the charge 
fractionization of type B. We pointed out that the topology that has been mentioned in 
the literature in connection with type B is of a trivial nature. 

The vacuum charge, which is the difference between the charge of the perturbed 
vacuum and that of the free vacuum, is in the form of the difference between two 
divergent series. The result depends on how the difference is calculated. We proposed 
a natural way of resolving this ambiguity. For the negative energy states of the free 
vacuum, we assume a Fermi energy Ep < 0. All the energy levels between and Ep 
are occupied and those below Ep are empty. When the external potential is switched 
on the Fermi energy shifts, Ep —>■ E' F but the one-to-one correspondence between the 
Es and E's is maintained. We eventually let Ep — > — oo and E' F — > — oo. This device is 
instrumental in clearly keeping track of the number of the particles in the vacuum. 

In Sec. Ill we examined a model that extends to the entire one-dimensional space. 
We divided the total charge of the vacuum Q into two parts, Q c and Qd- The Q c is 
the part of the charge that has been examined extensively in the literature. In order 
to see the charge conservation Q = Q c + Qd = 0, it is crucial to include Qd that is 
related to the change in the Fermi energy due to the external potential. Each of Q c and 
Qd may take a fractional value but we pointed out that this is a mathematical artifact 
that arises from replacing J2e with (R/tt) J dk. In this continuum limit the discrete 
energy levels become effectively smeared out. We emphasized that the level spacing and 
the energy shift in this situation are of the same order of magnitude. The number of 
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smeared out energy levels in a specified range becomes fuzzy. This artifact, however, 
does not jeopardize the conservation of the total charge. 

The model that we examined in Sec. Ill can be taken as a one- dimensional 
simulation of the relativistic shell model of the atomic nucleus. If we take the fractional 
charge literally, this leads to bizarre consequences as we discussed towards the end of 
Sec. III. If the shell model potential is constructed in a self-consistent manner, we 
expect that the diffuse distribution Pd{x) becomes confined to a finite region around the 
nucleus. The total nuclear charge is not fractionalized. 

In Sec. IV we examined the vacuum charge of a finite bag confined within \x\ < R. 
Since R is finite and not very large, we do not replace Y^e with / dk. Nevertheless a 
fractional charge can appear when we use the summation method, Eq. (14. 8p . that is 
related to the split-point regularization of the current operator. This occurs in the MIT 
bag model as well as in the chiral bag model. We concur with Banerjee [25] in arguing 
that this outcome of Eq. (14. 8 P is unnatural. From our point of view, no fractionization 
of the charge or the quark number appears in the bag models. 

What we have discussed in this paper may have relevance to the charge 
renormalization of quantum electrodynamics (QED). Suppose a test charge eo is placed 
in the vacuum. The eo is renormalized to e, 



The Ward identity leads to Z\ = Z 2 . We are still left with Z 3 which is due to the 
vacuum polarization; see, e.g. Ref. [26]. If we calculate Z 3 by perturbation theory, it is 
represented by a divergent integral. If we cut-off the integral, Z 3 depends on the cut-off 
parameter. Depending on the choice of the cut-off parameter, Z 3 can take any frac- 
tional value. This was probably the first example in which a fractional charge appeared. 
The process that is responsible for the charge renormalization is typically through the 
electron-positron pair creation in the virtual state. Because the charge of this pair is 
zero, however, it cannot change the charge of the system unless part of the charge some- 
how disappears. The mechanism of Z 3 ^ 1 is not clear. If we consider the same problem 
but assuming that the space is confined in a finite cavity, the situation will be similar to 
the one that we discussed in Sec. IV. Whatever happens to the system, the total charge 
is confined within the cavity and hence it is conserved. This seems to mean Z A — \. If we 
let the radius of the cavity tend to infinity, we will still have Z 3 = 1. Would it be possi- 
ble to reformulate QED such that we do not have to consider the charge renormalization? 

This work was supported by the Natural Sciences and Engineering Research Council 
of Canada. 
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